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1.0  ABSTRACT 


This  report  is  part  of  a  research  project  conducted  at  the 
Woods  Hole  Oceanographic  Institution  to  improve  the  flight  characteristics 
of  CTD*  instrument  packages.  Improvement  of  these  cable  lowered 
instrument  packages  could  allow  their  use  in  more  severe  weather  conditions. 
It  could  improve  the  quality  of  the  measurements. 

This  report  presents  the  development  of  a  simplified  mathematical 
model  of  the  CTD  package  flight  characteristics.  This  computer  model 
was  exercised  to  perform  a  sensitivity  analysis  of  different  versions 
of  CTD  packages. 

Part  of  the  research  project  includes  scale  model  testing.  The 
second  part  of  the  report  discusses  pertinent  flow  similarity  criteria 
and  proposes  a  scheme  for  building  a  CTD  half  scale  model. 

Finally,  recommendations  to  improve  the  hydrodynamic  behaviour  of 
the  present  CTD  configuration  are  summarized  at  the  end  of  the  report. 


*CTD  stands  for  Conductivity,  Temperature  and  Depth 


2.0  BACKGROUND 


The  main  users  of  CTD  instruments  are  physical  oceanographers 
who  wish  to  obtain  density  gradients  and  acousticians  who  wish  to  obtain 
sound  velocity  profiles.  The  CTD  measures  conductivity,  temperature  and 
pressure  from  which  the  density  and  sound  velocity  can  be  determined. 

A  typical  CTD  package  contains  three  instruments;  CTD,  rosette  and  a 
pinger,  all  mounted  within  a  protective  cage.  The  CTD  measures  the 
conductivity  (to  derive  salinity) ,  temperature  and  depth  of  the  package 
as  it  is  lowered  by  an  armored  electrical  cable.  The  rosette  consists  of 
either  twelve  or  twenty-four  water  sampling  bottles,  each  with  a  capacity 
of  1.2  to  4.0  liters.  The  acoustic  pinger  provides  a  measurement  of  the 
distance  off  the  bottom  in  order  to  obtain  samples  close  to  the  bottom  and 
to  prevent  the  CTD  package  from  impacting  the  bottom. 

Knowledge  of  the  density  structure  obtained  from  measurements  of 
salinity  and  temperature  at  a  specific  location  allows  the  oceanographer 
to  calculate  the  large-scale  circulation  of  water  masses  in  the  ocean. 

Water  samples  permit  a  check  on  the  accuracy  of  the  CTD  measurements. 

Samples  are  also  used  to  study  the  biological  and  chemical  properties  of 
the  ocean.  The  water  must  flow  freely  through  the  water  bottles  in  order 
to  prevent  any  entrapment  of  water  as  the  package  is  lowered  to  the  bottom. 

Only  modifications  which  preserve  or  improve  the  accuracy  of  the  measure¬ 
ment  described  above  should  be  implemented. 

The  following  scenario  describes  a  typical  CTD  cast.  The  research 
vessel  arrives  at  the  location  where  the  cast  is  to  be  made.  The  instrument 
is  eased  over  the  side  as  the  ship  heaves  to.  With  present  configurations, 
the  CTD  is  lowered  at  70  meters/minute  in  calm  seas  and  at  slower  speeds  in 
rougher  seas.  The  CTD  is  lowered  to  just  above  the  bottom  and  then  raised. 

Water  samples  are  collected  on  the  way  to  the  surface  at  predetermined 
depths.  To  take  a  water  sample,  the  CTD  is  brought  to  the  desired  depth 
and  allowed  to  equilibrate  with  the  surrounding  water  for  five  minutes. 

A  signal  is  then  sent  from  the  ship  causing  one  bottle  to  close.  The  CTD 
is  then  raised  to  the  next  sample  point  and  the  sequence  continued  until  the 
CTD  package  reaches  the  surface.  The  package  is  then  placed  back  on 
board,  the  water  samples  transferred  and  the  instruments  readied  as  the 
ship  steams  to  the  next  cast  site. 

As  can  be  inferred  from  the  above  scenario,  a  significant  length 
of  time  is  required  for  one  cast  to  be  made,  six  to  eight  hours  being 
typical.  With  the  high  cost  of  ship  time,  modifications  which  could  reduce 
the  amount  of  time  per  cast  would  be  beneficial.  The  primary  reason  for 
limiting  instrument  lowering  speed  is  to  insure  that  the  cable  supporting 
the  CTD  never  goes  slack,  which  causes  cable  kinking  and  electrical  failure, 

(Reference  2).  Increase  in  the  terminal  velocity  would  permit  a  faster 
lowering  speed  and  thereby  reduce  the  time  necessary  for  each  cast. 

One  other  characteristic  of  CTD  flight  concerns  the  stability  of 
the  motion.  It  is  possible  that  CTD's  in  their  present  configurations  go 
through  a  kiting  or  tumbling  motion  as  they  are  lowered.  The  effects  of 
such  motion  on  the  quality  of  scientific  data  is  not  presently  known.  One 
goal  of  this  research  project  is  to  investigate  the  stability  of  CTD  packages 
and  to  suggest  modifications  to  improve  the  flight  characteristics. 
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3.0  STATIC  ANALYSIS 

3 . 1  Static  Stability 

Static  stability  concerns  the  response  of  a  CTD  package 
when  it  is  disturbed  from  its  steady  state  vertical  flight.  The  CTD 
package  is  statically  stable  if  it  has  a  tendency  to  return  to  its  steady 
state  path  of  flight  after  removal  of  the  perturbation. 

Dynamic  stability  involves  the  path  the  perturbed  CTD  follows  after 
introducing  the  perturbation.  The  CTD  is  dynamically  stable  if  the 
oscillations  the  CTD  goes  through  after  perturbation  have  a  decaying 
amplitude.  A  short  discussion  of  CTD  dynamic  stability  is  presented  in 
Section  4. 

Ideally  a  good  CTD  instrument  package  should  be  both  statically  and 
dynamically  stable. 

3 . 2  CTD  Geometry 

CTD  instrument  packages  used  by  the  oceanographic  com¬ 
munity  differ  widely  in  both  the  number  of  components  used  to  make  up  the 
package  and  in  the  specifics  of  each  individual  component. 

Components  commonly  encountered  in  CTD  packages  used  at  the  Woods 
Hole  Oceanographic  Institution  include: 

1.  Neil  Brown  Instrument  Systems  CTD  with  24  inch  pressure  case. 

2.  General  Oceanics  (24) -1.2  liter  water  bottle  rosette. 

3.  WHOI  designed  3/4  inch  diameter  stainless  steel  protective  frame. 
Typically  the  frame  diameter  is  32  inches  and  its  height  is 

73  inches.  The  CTD  is  mounted  below  the  rosette  and  both  instruments 
are  surrounded  by  the  frame  as  can  be  seen  in  Figure  1. 

4.  An  optional  acoustic  pinger,  such  as  Benthos  type  2216,  which 
can  be  mounted  next  to  the  CTD  instrument  or  strapped  on  the 
outside  of  the  protective  frame. 

Mechanical  properties  of  actual  components  such  as  weight,  buoyancy, 
centers  of  gravity  and  buoyancy  required  for  this  analysis  have  been 
established  either  by  measurements  and/or  computation. 

3.3  Analytical  Model 
3.3.1  Assumptions 

Many  assumptions  were  necessary  for  developing  a 
simplified  mathematical  model  of  the  CTD  package.  First  the  complex  geometry 
of  the  package  had  to  be  simplified.  Components  were  assumed  to  be  of  simple 
shape  with  predictable  hydrodynamic  characteristics.  For  example  the  rosette 
has  been  modeled  as  a  blunt  cylinder  23  inches  in  diameter  and  30  inches  long. 
Figure  2  shows  the  simplified  geometry  of  the  package  main  components. 

Next  certain  assumptions  had  to  be  made  regarding  drag  coefficients. 

In  the  absence  of  experimental  values,  the  selection  of  drag  coefficients 
was  based  on  the  simplified  shape  of  each  component  for  both  laminar  and 
turbulent  flow  conditions.  For  the  CTD,  rosette,  and  pinger,  the  flow  based 
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on  Reynolds  number  was  just  below  the  transition  region  for  smooth  bodies. 
Turbulent  values  of  Cq  were  used,  however,  to  account  for  the  many  flow 
disturbing  appendages  of  these  components.  For  the  CTD  frame,  the  Reynolds 
number  was  much  below  transition  and  both  laminar  and  turbulent  values  of 
Cp  were  used. 

With  these  assumptions  cross  section  areas  and  drag  coefficients  for 
normal  and  tangential  flow  could  be  determined  for  the  different  CTD 
components.  Table  1  "Drag  Coefficients  and  Frontal  Areas  of  CTD  Components" 
presents  a  summary  of  these  values. 

3.3.2  Forces  at  Play 

The  forces  which  act  on  the  CTD  package  as  it  is 
steadily  lowered  in  the  ocean  are:  gravity,  buoyancy,  hydrodynamic  drag  and 
cable  tension. 

The  gravity  and  buoyancy  forces  are  easy  to  calculate  from  the  geometry 
and  weight  distribution  of  the  CTD  package  components. 

Drag  forces  are  more  difficult  to  predict.  They  are  resolved  into 
normal  and  tangential  components,  using  areas  and  drag  coefficients  of  the 
CTD  simplified  model  previously  described. 

The  tension  in  the  cable  supporting  the  CTD  package  can  be  inferred 
from  the  cable  angle,  the  package  weight  and  inclination,  and  the  drag 
forces  for  the  prevailing  lowering  speed. 

3.3.3  Mathematical  Model 


The  steps  that  the  mathematical  model  follows  to 
investigate  the  static  stability  of  a  particular  CTD  package  configuration 
include: 


o  Assume  an  initial  CTD  package  pitch  angle, 
o  Assume  a  lowering  speed. 

o  Resolve  the  resulting  drag  forces  into  normal  and  tangential 
components . 

o  Compute  the  moment  with  respect  to  the  CTD  package  center  of 
gravity (c.g. )due  to  the  hydrodynamic  drag  forces. 

o  Assume  a  cable  angle 

Derive  the  cable  tension.  The  formula  used  is: 


o 


DRAG  COEFFICIENTS  AND  FRONTAL  AREAS  OF  CTD  COMPONENTS 

TABLE  1 


I  1 .  TANGENTIAL  FLCW 
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where 


T 


B 

DTi 

DNi 

9 


Cable  tension 

Angle  between  cable  and  the  vertical 

Package  weight 
Package  Buoyancy 

y  I  f 

Tangential  drag  of  component  t- 
Normal  drag  of  component  " L* 

Package  inclination  angle  (pitch) 


o  Compute  the  moment  with  respect  to  the  CTD  t.g.  due  to  the  tension 
force. 

o  Compute  the  moment  with  respect  to  the  CTD  c.g.  due  to  the  package 
buoyancy  forces. 

o  Sum  the  moments. 

By  convention,  if  the  sum  is  negative  the  overall  moment  is  a  righting 
moment.  The  package  will  have  a  tendency  to  return  to  the  vertical  thus 
decreasing  its  pitch  angle.  The  configuration  can  be  considered  statically 
stable. 

If  the  sum  is  positive  the  overall  moment  is  a  capsizing  moment.  The 
pitch  angle  will  have  a  tendency  to  increase.  The  configuration  is  no 
longer  statically  stable. 

Figure  3  illustrates  how  a  typical  CTD  package  is  modeled.  Corresponding 
forces  and  sign  conventions  are  depicted  in  Figure  4. 

A  computer  program  called  CTDFL1  has  been  written  to  calculate  the 
resultant  moment  for  various  lowering  speeds.  Program  inputs  include  the 
physical  and  geometrical  characteristics  of  the  CTD  package  components, 
the  angle  between  the  lowering  cable  and  the  vertical,  and  the  package 
inclination.  Program  outputs  include  cable  tension,  tangential  and  normal 
drag  forces  and  overall  resultant  moment.  The  lowering  speed  is  increased 
until  the  tension  in  the  cable  vanishes.  The  package  is  then  assumed  to 
have  reached  terminal  velocity.  The  speed  at  terminal  velocity  is  also  a 
program  output.  Appendix  A  is  a  description  of  CTDFLI. 

3.4  Sensitivity  Analysis 

The  program  CTDFLI  was  exercised  to  study  the  static 
stability  of  typical  and  modified  CTD  instrument  packages. 

The  stability  of  typical  CTD  packages  was  first  established.  Changes 
were  then  made  to  the  base  line  configurations  and  the  resulting  effects 
were  investigated.  Changes  included  increased  weight,  reduced  drag,  and 
improved  package  symmetry.  The  sensitivity  of  the  model  to  changes  in 
tension  and  tension  angle  was  also  investigated  using  the  standard  CTD  package. 
Finally  several  runs  using  combinations  of  the  changes  mentioned  were  made  in 
an  attempt  at  finding  optimal  configurations.  Table  2  lists  all  the  study 
cases  investigated. 


CG  -  Center  of  Gravity  0  - 

CB  -  Center  of  Buoyancy  P  - 

Vn  -  Normal  Drag  Force  F  - 

Dj  -  Tangential  Drag  Force  R  - 

T  -  Tension  in  cable  CTD- 


Overall 
Pinger 
Frame 
Rosette 
CTO  Case 


Figure:  3 
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configuration 
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Same  as  A.O,  cable  angle  *  10 


Same  as  A.O,  cable  angle  =  20 


Same  as  A.O,  cable  angle  =  30 


Same  as  C.0,  cable  angle  -  10 


Same  as  C.0,  cable  angle  ■  20 


Same  as  C.0  cable  angle  =  30 
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steel. 


A.0.W2  Same  as  A.O,  add  50  lbs  to  CTD 


A.0.W3  Add  one  inch  lead  ring  attached  to  frame  lower  ri 


A.0.W4  Add  two  inch  lead  ring  attached  to  frame  lower  rin 


odif ications 


C.0.W3  Same  as  C.0,  except  add  1  in.  lead  ri 


C.0.W4  I  Same  as  C.0,  except  add  2  in.  lead  ri 


C.0.D1  Same  as  C.0,  except  frame  drag  reduced 


C.0.D2  Same  as  C.0,  except  frame  drag  reduced  to 


C.0.D3  Same  as  C.0,  except  frame  drag  reduced  to  1/8 


C.0.D4  I  Same  as  C.0,  except  frame  size  reduced 


n  C.0 


[configuration 


F.O  Same  as  D.0,  except  drag  reduced  %  and 


1  in.  lead  ring  attached 


timum 


Same  as  above  except  small  WHOI  frame 
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A  description  of  the  configurations  studied  and  a  review  of  the  results 
of  the  computer  analysis  are  hereafter  presented. 

3.4.1  Present  Configurations 

Three  instrument  packages,  typical  of  those 
commonly  encountered  in  practice  were  first  modeled.  These  configurations 
were  as  follows: 

o  Configuration  A.  Neil  Brown  24  in.  CTD,  General  Oceanics  24-1.2  liter 
rosette,  WHOI  standard  protective  frame. 

o  Configuration  B.  Same  as  A  except  for  Benthos  pinger  strapped 
to  CTD. 

o  Configuration  C.  Same  as  A  except  for  Benthos  pinger  strapped 
to  frame. 

These  three  basic  CTD  package  configurations  were  used  to  determine  baseline 
data  for  the  sensitivity  analysis  and  to  verify  the  validity  of  the  model. 

The  terminal  velocity  obtained  for  these  three  configurations  ranged  from 
6.26  ft/sec  to  6.95  ft/sec.  An  actual  drop  test  conducted  in  July  1979  with 
a  configuration  similar  to  B  yielded  a  terminal  velocity  of  6.85  ft/sec. 

Good  agreement  seems  therefore  to  exist  between  the  computed  and  the  measured 
terminal  velocities. 

Righting  and/or  capsizing  moments  were  computed  for  the  three  configurations 
assuming  their  pitch  angle  to  vary  from  0  to  45  degrees.  A  typical  curve 
of  moments  as  a  function  of  package  speed  is  presented  in  Figure  5. 

To  better  understand  this  curve  let  us  first  consider  what  happens 
to  configuration  A  as  it  changes  inclination  (pitch)  and  falling  speed. 

The  external  forces  acting  on  the  package  inclined  at  some  angle  theta 
are  shown  in  Figure  6.  Some  of  these  forces,  drag  and  tension  for  example, 
vary  as  a  function  of  speed  and  their  relative  contribution  to  the  overall 
moment  needs  some  reflection. 

Since  configuration  A .0  is  axisymmetric  no  moment  can  exist  at  zero 
pitch  angle,  thus  the  zero  moment  value  shown  for  the  first  moment  curve. 

The  first  point  on  the  second  curve  gives  the  value  of  righting  moment  at 
speed  zero  and  five  degrees  of  pitch  angle.  This  static  moment  is  simply 
the  sum  of  two  moments.  The  first  moment  is  due  to  the  normal  component  of 
the  tension  force.  For  this  configuration,  the  tension  acts  vertically  so 
at  5°  of  pitch  this  force  creates  a  righting  moment  about  the  CG.  Both  these 
contributions  increase  as  the  pitch  angle  increases. 

As  the  package  acquires  speed,  drag  force  comes  into  play.  Normal 
and  tangential  drag  forces  are  computed  for  each  CTD  component  based  on  its 
drag  characteristics.  These  forces  are  then  applied  at  the  center  of 
gravity  of  each  of  the  CTD  components.  For  configuration  A  there  are  three 
components  as  described  previously.  If  only  drag  forces  were  present,  the 
righting  moment  curve  for  each  pitch  angle  would  simple  be  a  parabola  (since 
drag  «c  v*  )  with  a  zero  velocity  value  of  zero.  Adding  the  buoyancy  force 
in  the  model  just  shifts  these  curves  down  by  an  amount  equal  to  the  static 
moment  described  above  without  tension.  It  is  advisable  to  keep  this  parabolic 
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T  -  Tension 
Dfi  -  Normal  Drag 
Dj  -  Tangential  Drag 
B  =  Buoyancy 
W  =  Weight 

R  =  Rosette 
F  =  Frame 

CTD  s  CTD  Instrument 
CGq  s  Overall  C.G. 


Figure:  6 


curve  in  mind  when  looking  at  the  additional  effect  the  cable  tension  has 
on  stability. 


As  described  earlier  in  the  report  the  tension  in  the  cable  is 
determined  by  taking  the  difference  between  the  submerged  weight  of  the 
package  and  the  vertical  components  of  drag.  At  zero  speed  the  magnitude 
of  the  tension  force  equals  the  submerged  weight  of  the  package.  As  the 
speed  increases  the  drag  forces  increase  with  the  velocity  squared  so  that 
the  tension  decreases  in  a  corresponding  way.  At  terminal  velocity  the 
tension  is  zero  and  the  righting  moment  reaches  its  smallest  value.  For 
all  pitch  angles  greater  than  zero,  the  righting  moment  is  a  maximum  at 
zero  speed  and  decreases  in  a  parabolic  fashion  as  the  speed  increases  (due 
to  decreasing  tension).  Once  terminal  velocity  is  reached  tension  has  no 
more  effect.  From  the  graph  it  is  readily  seen  that  at  speeds  below 
terminal  velocity  the  tension  force  has  the  greatest  effect  on  the  righting 
moment . 

Figures  7  and  8  show  the  stability  curves  for  the  next  two  basic 
configurations  studied,  B.O  and  C.O.  Both  these  configurations  are  not 
axisymmetrical,  configuration  C.O  with  a  pinger  strapped  on  the  outside 
of  the  frame  being  more  asymmetrical  than  B.O. 

It  can  be  seen  from  these  graphs  that  the  static  stability  varies 
markedly  with  symmetry.  Configuration  A.O  is  stable  throughout  a  range  of 
pitch  angles  of  0°  to  45° ,  indicating  that  the  CTD  when  perturbed  would 
tend  to  right  itself  to  zero  pitch  angle.  On  the  other  hand  configurations 
B.O  and  C.O  are  both  unstable  at  zero  pitch  angle  with  C.O  having  the 
largest  overturning  moment  as  expected.  The  consequences  of  these  instabilities 
in  actual  flight  are  unknown  at  present  but  several  things  are  obvious. 
Configurations  which  are  stable,  fall  with  zero  pitch  angle,  should  have 
a  higher  terminal  velocity  since  the  projected  area  of  the  package  is  a 
minimum.  Packages  with  overturning  moments  may  tend  to  oscillate  throughout 
the  range  of  instability  causing  kiting  or  tumbling. 

As  part  of  the  baseline  case  studies,  the  sensitivity  of  stability 
to  cable  angle  was  next  investigated.  The  symmetrical  A.O  and  the 
asymmetrical  C.O  configurations  where  subjected  to  cable  pull  angles  of 
10,  20  and  30  degrees  respectively.  Figures  9  to  14  show  the  stability 
curves  thus  obtained.  These  results  indicate  that  increasing  the  pull  angle 
causes  a  reduction  in  the  righting  moment,  and  thus  enhances  the  potential 
of  oscillations  of  a  lowered  CTD  package.  It  is  not  known  if  fluctuations  in 
tension  angles  occur  at  sea  although  it  has  been  hypothesized  (Reference  2). 

3.4.2  Symmetry  Modification 

To  test  the  sensitivity  of  the  model  to  symmetry, 
configurations  B.O  and  C.O  were  modified  to  make  them  axisymmetric . 
Configurations  D.O  and  E.O  were  obtained  by  simply  adding  another  pinger  - 
or  pinger  casing  -  opposite  the  original  one.  The  effect  on  stability  was 
immediate,  both  configurations  were  stable  throughout  the  range  of  pitch 
angles  previously  tested.  Given  the  mathematical  model  used  in  CTDFCI, 
this  result  was  to  be  expected.  The  stability  curves  of  D.O  and  E.O  are 
shown  in  Figures  15  and  16. 
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The  behavior  of  a  real  CTD  package  in  the  ocean,  however,  may  not 
have  the  same  response  although  it  seems  logical  that  a  symmetric  body 
would  fall  vertically.  An  added  bonus  to  the  flight  characteristics  of 
both  configuration  D  and  F  was  an  increase  in  terminal  velocity  of  .25  ft/sec. 
Modifying  a  package  to  make  it  synmetric  should  improve  the  stability 
greatly  and  it  is  recommended  that  in  the  future  more  concern  be  placed  with 
making  CTD  packages  as  syrmnetric  as  possible. 

3.4.3  Weight  Modifications 

The  sensitivity  of  the  model  to  the  judicious 
addition  of  weight  was  studied  next  using  again  both  the  A (symmetrical) 
and  the  C (asymmetrical)  configurations.  Weight  changes  were  introduced  as 
follows:  In  configuration  A.0.W1  the  bottom  ring  of  the  standard  WHOI  frame 
is  made  of  solid  steel  instead  of  tubing.  In  A.0.W2  50  lbs  are  arbitrarily 
added  to  the  CTD  instrument.  In  A.0.W3,  W4  and  C.0.W3,  W4  lead  rings  of  one 
and  two  inches  in  diameter  are  fastened  to  the  frame  lower  ring.  These 
changes  could  be  easily  implemented  in  practice.  The  increase  in  weight 
and  terminal  velocity  that  these  changes  caused  are  hereafter  tabulated. 


TABLE  3 


Weight  Modifications 


Case  # 

Original 

Weight 

Added 

Weight 

Terminal  Velocity  I 

Ft/Sec 

%  Increase 

A.0.W1 

400 

15 

6.69 

3.8 

A.0.W2 

400 

50 

7.16 

11.2 

A.0.W3 

400 

34.5 

6.91 

7.3 

A.0.W4 

400 

125.0 

7.94 

23.3 

C.0.W3 

455 

34.5 

7.38 

6.3 

C.0.W4 

455 

125.0 

8.31 

19.7 

Stability  curves  for  these  six  configurations  appear  on  Figures  17  to  22. 

The  major  effect  of  weight  addition  appears  to  be  a  substantial  increase  in 
package  terminal  velocity.  Stability  improvements  appear  less  obvious. 

3.4.4  Drag  Modifications 

At  present,  CTD  packages  are  anything  but  stream¬ 
lined.  Improvement  in  the  drag  characteristics  of  a  CTD  configuration  by 
adding  fairings,  shrouds,  varying  dimensions,  etc.,  should  improve  the 
terminal  velocity  substantially.  To  determine  the  sensitivity  of  the  model- 
to-drag  modifications  four  runs  were  made.  Three  of  these  investigate  the 
improvement  in  terminal  velocity  due  to  reducing  the  drag  on  the  WHOI  CTD 
frame  by  1/2,  1/4  and  1/8.  At  this  stage,  the  exact  mechanism  for  reducing 
the  drag  to  these  levels  has  not  been  determined  although  fairings  would  be 
a  good  candidate.  The  fourth  modification  involved  substituting  the  standard 
WHOI  CTD  frame  with  one  made  of  3/8  in.  pipe  instead  of  3/4  in.  pipe.  The  increase 
in  terminal  velocities  introduced  by  these  changes  are  hereafter  tabulated. 
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TABLE  4 

Drag  Modifications 


Case  # 

Drag  Reduction 

Terminal 

Ft/Sec 

Velocity 
°L  Increase 

C.0.D1 

1/2 

7.94 

14.4 

C.0.D2 

1/4 

8.63 

24.4 

C.0.D3 

1/8 

9.06 

30.5 

C.0.D4 

Small  Frame 

8.44 

21.6 

Figure  23  shows  the  stability  curves  for  zero  pitch  angle  of  these  four  cases. 

These  results  clearly  indicate  that  drag  reduction,  even  as  little  as 
greatly  improves  the  package  terminal  velocity.  A  judicious  combination  of 
weight  increase  and  drag  reduction  can  only  benefit  the  performance  of  future 
packages. 


3.4.5  Optimal  Configurations 

Based  on  the  results  of  the  previous  four  sections, 
two  runs  (F.O  and  G.O)  were  made  in  an  attempt  to  optimize  both  stability  and 
terminal  velocity.  Both  configurations  use  configuration  B.O  as  a  basis 
with  the  addition  of  another  pinger  to  make  it  symmetric,  a  one  inch  lead 
ring  to  add  weight,  and  a  reduction  in  drag  of  the  WHOI  frame  by  Con¬ 

figuration  G.O  differs  from  F.O  described  above  in  that  the  small  diameter 
WHOI  frame  is  used.  Stability  curves  for  these  two  last  runs  are  shown 
in  Figures  24  and  25.  These  curves  show  a  pronounced  increase  in  stability. 
The  improvement  in  terminal  velocity,  shown  in  the  table  hereafter,  is  even 
more  impressive. 


TABLE  5 

"Optimal"  Cases 


Case  # 

Terminal  Velocity 
Ft/Sec 

°L  Increase 

B.O 

6.94 

0 

G.O 

9.44 

36.0 

F.O 

8.84 

27.4 

These  values  should  not  be  thought  of  as  final,  however,  as  further  study 
should  investigate  ways  in  improving  these  numbers  even  more. 
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4.0  DYNAMIC  ANALYSIS 


tamic  Stability 


The  dynamic  stability  of  an  object  moving  at  a  constant 
velocity  involves  the  path  an  object  takes  around  the  steady  motion 
equilibrium  position  after  a  perturbation.  At  the  outset  of  this  study, 
it  was  thought  that  a  dynamic  stability  analysis  would  be  required  to 
determine  the  flight  characteristics  of  CTD  packages.  A  literature  survey 
was  conducted  and  the  preliminary  equations  of  dynamic  stability  were 
derived,  (Appendix  B) . 


To  conduct  a  meaningful  dynamic  stability  analysis,  the  values  of 
the  stability  derivatives  are  required.  Stability  derivatives  are  partial 
differentials  which  express  how  a  force  or  moment  changes  with  a  differen¬ 
tial  change  in  position,  velocity  or  acceleration  from  steady  motion.  A 
CTD  package  is  a  complicated  body  of  several  components  and  as  such  there 
is  no  analytical  method  to  determine  the  true  values  of  the  stability 
derivatives.  Values  of  stability  derivatives  can  be  inferred  if  enough 
simplifications  are  made  but  the  hope  that  the  resultant  model  replicates 
reality  is  remote. 


The  other  method  to  determine  the  values  of  the  stability  derivatives 
would  be  to  conduct  a  model  test  in  a  wave  tank  or  water  tunnel.  The 
CTD  model  would  have  to  be  fully  outfitted  to  measure  the  forces  and 
moments  in  the  principal  directions  and  to  be  able  to  be  perturbed  from 
its  steady  motion.  Such  a  test  program  is  beyond  the  scope  of  this  study. 


Vortex  Shedding 


It  has  been  inferred  that  alternate  vortex  shedding  on 
either  side  of  a  horizontally  mounted  CTD  (or  nephelometer)  may  contribute 
to  an  instability  in  CTD  flight  and  help  cause  large  scale  oscillations 
or  kiting*.  For  a  smooth  cylinder  vortex  shedding  would  definitely  be 
occurring.  The  effect  of  this  small  alternating  hydrodynamic  force  should 
not  be  significant  unless  of  course  its  frequency  is  close  to  the  pendulum 
frequency  of  the  CTD  package.  The  natural  frequency  of  oscillations  of 
the  CTD  package  about  its  point  of  attachment  to  the  lowering  cable  has  been 
estimated  to  be  .54  (Hz)  in  air.  The  shedding  frequency  of  the  horizontally 
mounted  CTD  (7"  diameter)  is  found  from 


where  V  is  the  lowering  speed  (ft/sec) 
and  D  is  the  CTD  diameter  (ft). 

♦Footnote.  Verbal  communication,  R,  Reiniger,  Bedford  Institute  of  Oceanography, 
Halifax,  Nova  Scotia,  CANADA. 
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At  a  lowering  speed  of  1.5  ft/sec  (27  meters/minute)  the  shedding 
frequency  would  be  close  to  the  CTD  pendulum  frequency  and  resonance  would 
occur  with  possible  side  motion  dependent  on  the  damping  in  the  system. 
However  most  lowerings  are  done  at  higher  speeds,  typically  60  to  70 
meters/minute.  At  70  m/tnin  the  frequency  of  the  vortex  shedding  induced 
by  the  horizontally  mounted  CTD  is  found  to  be  1.3  Hz.  As  can  be  seen 
from  the  simple  calculations  presented  in  Appendix  C,  the  package 
oscillations  resulting  from  this  vortex  shedding  would  have  a  very  small 
amplitude  (less  than  a  degree).  Furthermore  in  actual  practice  the  clamps 
which  hold  the  horizontally  mounted  instruments  cause  turbulence  which 
should  reduce  the  likelihood  of  periodic  vortex  shedding  substantially. 

It  thus  appears  that  vortex  shedding  is  not  a  serious  problem,  if  it  exists 
at  all. 

5.0  MODEL  TESTING 


To  actually  determine  the  flight  characteristics  of  existing 
and  modified  CTD  packages  two  series  of  tests  will  be  conducted  as  part 
of  the  current  CTD  hydrodynamics  study. 

Model  testing  will  be  first  conducted  in  the  50  foot  wide  and  100 
foot  deep  tank  of  the  Naval  Surface  Weaoons  Center  in  Silver  Spring, 

Maryland.  Full  scale  tests  will  be  later  conducted  at  sea,  using  an 
actual  CTD  package  specially  instrumented  to  monitor  the  cable  tension 
at  the  package  and  the  package  flight  pattern. 

The  objectives  of  the  scale  model  tests  will  be  to  observe  and 
record  the  flight  stability  of  different  scaled  down  versions  of  CTD  packages 
as  they  are  cable  lowered  or  free  fall  to  the  bottom  of  the  tank.  The 
rationale  for  the  selection  of  the  scaling  factors  and  the  resulting 
geometry  and  materials  for  model  fabrication  are  hereafter  explained. 

5 . 1  Dimensional  Analysis 

The  feasibility  for  constructing  a  smaller  scale  model 
of  the  CTD  package  is  based  on  dimensional  analysis.  The  use  of  a  scale 
model  is  desired  since  the  model  would  be  smaller,  easier  to  handle,  less 
expensive  to  repair  if  damaged  and  easily  modified.  The  liability  of  a 
model,  however,  is  the  uncertainty  involved  over  whether  the  model 
duplicates  the  flight  of  the  full-scale  CTD.  Proper  hydrodynamic  scaling 
and  geometric  similarity  between  model  and  full  scale  must  minimize  the 
scale  effects  between  the  two  packages. 

In  a  model  test  there  are  three  kinds  of  similarity  which  must  be 
satisfied: 

1.  Geometric  similarity. 

2.  Kinematic  similarity.  Streamline  pattern  in  the  model  must 
be  the  same  for  model  and  prototype. 

Dynamic  similarity.  Force  ratio  must  be  the  same  between 
model  and  prototype. 


3. 
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Geometric  similarity  is  relatively  easy  to  obtain.  It  reproduces 
not  only  the  shape  of  the  prototype  but  also  its  mass  and  displacement 
distribution  (center  of  gravity,  center  of  buoyancy,  moment  of  inertia, 
etc . ) . 

In  flow  regimes  where  the  drag  coefficients  are  highly  dependent 
on  the  Reynolds  number,  kinematic  similarity  will  dictate  that  the 
Reynolds  number  be  the  same  for  both  prototype  and  model.  Inasmuch  as 
the  Reynolds  number  ^  is  of  the  form 


where  K  is  the  inverse  of  the  fluid  kinematic  viscosity,  \/  is  the 
speed  of  the  object  and  J3  its  dimension,  an  equality  of  Reynolds 
numbers  for  model  and  prototype  tested  in  the  same  fluid  implies  that 


For  this  equality  to  be  maintained  the  speed  of  a  ^  scale  model  must 
therefore  be  twice  as  large  than  the  speed  of  prototype. 


Furthermore,  at  terminal  velocity  the  drag  of  a  free  falling  object 
must  equal  its  submerged  weight. 

The  familiar  equation  expressing  this  fact  is: 

5-  ±t>Q,s4Vln  W 


f 

c. 


where  ^  is  the  drag  force 

is  the  fluid  density 
is  the  drag  coefficient 
is  a  representative  area  of  the  object 

where ^  is  some  geometry  constant 
1/  is  the  terminal  velocity, 

and  w  is  the  immersed  weight,  that  is  the  difference  between 
the  object's  actual  weight  and  the  weight  of  the  water 
it  displaces. 


The  corresponding  equations  for  model  and  prototype  falling  at  terminal 
velocity  are: 


i  J c*» &  ^ 
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Now  for  geometric  and  kinematic  similarity  to  prevail. 


Cj>Ai  c 

3* 


and  therefore 

^  =  Wp 

One  may  note  that  if  then  the  condition  of  dynamic  similarity, 

namely 

Fz>+4  =r  *1>F> 

U/m  top 

is  automatically  satisfied. 

A  half  scale  model  which  falls  twice  as  fast  and  weighs  in  water  as  much 
as  the  prototype  is  not  very  desirable.  Fortunately  it  has  been  established 
that  the  drag  coefficient  of  blunt  objects  in  turbulent  flow  regimes  remains 
approximately  constant  over  a  large  range  of  Reynolds  numbers. 

This  being  the  case  and  provided  that  flow  turbulence  is  maintained  for 
the  model,  only  the  dynamic  similarity  will  be  used  as  our  second  criterion  - 
geometric  similarity  being  the  first  -  for  the  design  of  the  model. 

5.2  Half  Scale  Model  Design 

Maintaining  the  drag-to-submerged  weight-ratio  equal 
between  model  and  full-scale  requires  satisfaction  of  the  following  equations 

e  i~S>r> _ 

Wai  M/p 

or 

ton  Wj> 


Assuming  that  *  Cjp  ? 

and  that 

and  selecting  a  speed  of  model  fall  equal  to  the  speed  of  prototype  fall  - 
which  should  be  sufficient  to  maintain  turbulent  flow  conditions  -  yields 

k/fr  b  lA//> 

f 
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In  other  words,  a  %  scale  model  which  has  the  same  terminal  velocity  as 
the  full-scale  CTD  would  have  a  submerged  weight  equal  to  \  of  the  full- 
scale  value.  This  result  is  used  next  to  develop  specifications  for  a 
%  scale  model. 


The  geometry  of  the  scale  model  will  be  patterned  after  the  simplified 
geometry  used  in  the  mathematical  model.  The  different  CTD  components  will 
be  homogeneous  cylinders  with  one-half  the  length  and  one-half  the  diameter 
of  the  cylinders  they  represent. 


To  maintain  correct  submerged  weight  distribution  as  well  as  correct 
total  model  weight,  the  material  density  of  each  cylinder  must  be  calculated 
according  to  the  following  procedure: 


1. 

2. 


ft  •  * 


Establish  the  immersed  weight  of  the  actual  component  u  Up  CO 

//  <  • 

The  immersed  weight  of  the  model  component  C  is  then 

o)  =  wfcl)  /a 


3.  Compute  the  volume 


vau) 


p  •  ** 

of  the  model  component  t- 


The  density  of  the  model  component 


is  the 


the  density  of  water  (lbs/cu.ft). 


Results  obtained  with  this  simple  computation  procedure  are  sumnarized  in 
Table  6. 


TABLE  6 

Specifications  for  \  Scale  CTD  Model 


Component 

Modeled 

Cylinder  Dimensions 
(inches) 

(Diam  x  Length) 

Weight 

(lbs) 

Immersed 

Weight 

(lbs) 

Specific 

Weight 

(lbs/cu.ft. 

I 

CTD 

3.5  x  12.5 

29.26 

24.75 

415.63 

6 . 66 

Rosette 

11.5  x  15.0 

65.20 

17.50 

mam 

1.34 

Frame 

\  0  diam  tubing 

Ring  &  Verticals 
to  scale 

16.90 

13.50 

310.0 

4.97 

Pinger 

2.5  x  13 

12.24 

10.0 

349.71 

5.60 

Nephelo 

3.0  x  16 

22.87 

18.68 

—  ...  - 

349.71 

5.60 
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The  immersed  weight  of  a  CTD  prototype  package  made  of  a  CTD  instrument, 
a  rosette,  a  WHO  I  frame  and  one  pinger  is  263  lbs.  The  immersed  weight 
of  this  package  model  is  found  to  be  65.25  or  4.03  as  small  as  the 
immersed  weight  of  the  prototype. 

6.0  RECOMMENDATIONS  AND  CONCLUSIONS 

The  sensitivity  analysis  has  shown  that  substantial  improve¬ 
ment  in  CTD  performance  is  possible.  Symmetry  was  found  to  be  very 
important  to  the  stability  of  CTD  packages  and  this  result  should  be 
demonstrated  in  the  NSWC  tank.  The  reduction  of  drag  was  also  determined 
to  significantly  increase  the  terminal  velocity  of  CTD  configurations 
and  modifications  to  achieve  this  reduction  should  be  much  more  thoroughly 
studied  in  the  final  phase  of  the  project.  A  search  for  optimal  config¬ 
urations  combining  the  beneficial  effects  of  symmetry,  drag  reduction 
and  weight  addition  should  be  concentrated  on.  Scale  models  should  be 
constructed  to  use  during  the  test  phase  which  can  be  easily  modified. 

In  this  preliminary  study  a  simple  mathematical  model  was  developed 
and  used  to  study  the  hydrodynamic  characteristics  of  CTD  packages  and 
the  sensitivity  of  the  package  to  various  modifications.  The  form  of  the 
analytical  model  used  above  is  adequate  for  this  preliminary  study.  It  is 
hoped  that  a  much  more  detailed  analysis  be  conducted  to  more  accurately 
model  the  hydrodynamics  of  CTD  and  other  oceanographic  instrument  packages. 
A  fully  instrumented  model  with  which  to  obtain  stability  derivatives  would 
be  a  logical  first  step. 

Scale  model  tests  and  tests  performed  at  sea  with  a  full  scale 
model  should  confirm  the  trends  point  out  by  relatively  simple  and  more 
sophisticated  analyses. 

The  ultimate  reconmendation  will  be  to  sensibly  alter  the  existing 
configuration  and  to  design  a  new  package  which  would  greatly  improve 
the  present  performance. 
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APPENDIX  "A" 


"Description  of  Computer  Program  CTDFLI" 


Appendix  A 

Description  of  Computer  Probram  CTDFLI  7-28-80 

Name :  CTDFLI 

Type:  Main  Program 

Purpose:  Calculates  the  fluid  drag  forces  and  moments  acting  on  CTD 

packages  to  determine  the  stability  of  vertical  motion  as 
a  function  of  velocity  and  angle  of  orientation. 

Machine:  Xerox  Sigma  7 

Source  Language:  Xerox  Extended  FORTRAN  IV 

Program  Category:  Numerical  Model 
Description: 

The  program  computes  the  normal  and  tangential  drag  forces,  the 
pitching  moments  these  produce  and  the  overall  righting  moment  on  CTD 
configurations  as  a  function  of  vertical  velocity  and  angle  of  orientation. 

A  CTD  configuration  is  made  up  of  up  to  ten  different  components  of  known 
center  of  gravity  and  tangential  and  normal  drag  constants.  The  weight, 
buoyancy,  center  of  buoyancy,  cable  tension  angle  and  point  of  application 
of  the  cable  tension  must  also  be  known  for  each  CTD  configuration.  For 
each  specific  angle  of  orientation,  the  normal  and  tangential  drag  forces 
are  calculated  for  each  component  at  velocities  from  0  to  terminal  velocity 
at  0.5  ft/sec  increments.  The  pitching  moment  about  the  overall  centers  of 
gravity  of  a  configuration  due  to  the  drag  forces  is  then  calculated  as 
well  as  an  overall  righting  moment.  In  addition  to  the  pitching  moment, 
the  righting  moment  is  made  up  of  the  moments  caused  by  the  buoyancy  force 
and  by  the  cable  tension.  The  terminal  velocity  at  each  angle  of  orientation 
is  also  calculated  and  output. 

The  program  has  been  written  to  run  either  on-line  or  as  a  batch  job, 
with  the  primary  use  assumed  to  be  on-line.  The  user  must  supply  parameters 
defining  the  geometry  of  the  configuration  and  the  drag  characteristics 
of  each  component.  The  center  of  buoyancy  and  centers  of  gravity  of  each 
component  are  measured  relative  to  the  overall  center  of  gravity  of  each 
configuration.  The  overall  center  of  gravity  can  be  thought  of  as  the  origin 
of  a  body-fixed  coordinate  system  with  the  x-axis  positive  (tangential  direction) 
positive  downward  and  the  z-axis  (normal  direction)  positive  to  the  right. 

The  user  supplied  parameters  are: 

BUOYF  Buoyancy  force  of  overall  package,  lbs 

XCB  Distance  along  x-axis  to  center  of  buoyancy  (CB) ,  inches 

ZCB  Distance  along  z-axis  to  center  of  buoyancy,  inches 

WTAIR  Weight  in  air  of  CTD  package,  lbs 

PHI  Angle  tension  acts  measured  from  the  vertical,  degrees 

XTEN  Distance  along  x-axis  to  pt.  of  application  of  tension,  inches 

Distance  along  z-axis  to  pt.  of  application  of  tension,  inches 


ZTEN 
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For  each  component: 

XCGi  Distance  along  x-axis  to  Center  of  Gravity  (CG)  of  component  i, 

inches 

ZCGi  Distance  along  z=axis  to  CG  of  component  i,  inches 

CDANi  Normal  drag  constant  of  component  i 

CDATi  Tangential  drag  constant  of  component  i 

THETA  Pitch  angle;  angle  between  body-fixed  coordinates  and  the  vertical, 

degrees  (Up  to  8  different  pitch  angles  can  be  specified  per  run) 

When  all  the  parameters  have  been  input  the  calculation  sequence  begins 
for  the  pitch  angle  specified.  Forces  and  moments  are  calculated  for  31 
different  velocities  between  0  and  15  ft/sec  at  0.5  ft/sec  intervals.  The 
calculation  sequence  for  one  velocity  is  as  follows: 

1.  Break  velocity  into  normal  and  tangential  components 

VELN  =  VEL*  SIN (THETA) 

VELT  =  VEL*  COS (THETA) 

2.  For  each  component,  calculate  tangential  force: 

XI  =  CDAT(I)*  VELT*VELT 
normal  force: 

ZI  =  CDAN(I)  *  VELN*VELN 
and  pitching  moment: 

PMI  =  CDAN ( I ) *VE  LN*VE  LN*XCG ( I )  +  CDAT(I)*VELT*VELT*ZCG(1) 

3.  Sum  these  up  to  get  the  total  force  and  moments: 

XF  total  tangential  force,  lbs 

ZF  total  normal  force,  lbs 

PMDYN  pitching  moment,  ft/'lbs 

4.  Calculate  tension  in  the  cable: 

TENS0N(J)  =  1/C0S(PHI)*((SUBWT  -  XF  *  COS (THETA)  -  ZF*S IN (THETA ) ) 

5.  The  program  now  checks  to  see  if  terminal  velocity  has  been  reached. 
Terminal  velocity  occurs  when  the  tension  in  the  cable  is  zero.  If 
terminal  velocity  has  been  reached,  which  is  determined  when  the 
tension  goes  from  a  positive  to  a  negative  value  at  the  next  velocity 
increment,  the  terminal  velocity  is  calculated  by  interpolation  using 
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the  values  of  tension  at  the  present  and  previous  velocity  values. 

The  formula  used  is 

VTERM  =  TENSON(J-l) / (TENSON(J -1)  -  TENSON(J))  *  0.5  +  (VEL  -  .5) 

Once  terminal  velocity  has  been  reached  no  more  moment  calculations 
are  performed  for  this  pitch  angle. 

6.  Calculate  righting  moment: 

PMSTAT  =  PMDYN  +  BUOYF  *  SIN(THETA)  *  XCB  +  BUOYF*COS (THETA) *ZCB  + 

TENSON(J)*S IN (THETA-PHI )*XTEN  +  TENSON(J)*COS(THETA-PHI)*ZTEN 

Following  the  calculation  sequence,  the  results  can  be  printed  out. 

The  program  then  asks  if  another  run  using  a  different  pitch  angle 
is  wanted.  If  so,  a  new  value  of  pitch  angle  is  input  and  the 
calculation  sequence  is  repeated.  If  no  other  runs  are  needed, 
the  program  asks  the  user  if  plots  are  desired.  Plots  of  tangential 
drag  force,  normal  drag  force,  pitching  moment  and  righting  moment 
versus  velocity  can  be  obtained.  The  user  specifies  which  plot  he 
wants  drawn  and  the  title  of  the  plot.  Up  to  four  plots  can  be 
specified  per  run.  At  this  point,  the  program  stops.  The  details 
involved  in  having  the  plots  printed  are  described  in  the  Operating 
Instructions  below. 

Operating  Instructions 

I.  On-line  usage 

1.  To  start  the  program  respond  to  the  I  prompt  from  the  system 
by  entering: 

PLATE N_80 

This  coranand  allows  output  to  fit  on  paper.  The  computer  should 
respond  with  another  !  prompt,  then  type 

SET  F: 95/Filename 

This  command  creates  a  plotting  file  of  the  name  specified  which 
will  contain  the  plotting  information  the  user  specifies  while 
running  the  program.  The  computer  should  respond  with  another 
! prompt,  then  type. 

S_CTDRUN 

2.  Program  should  respond 

CTDFLI- -VERS ION  1.0  -  JULY  1980 
INPUT  BUOYANCY  FORCE  (LBS) 

? 
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Respond  by  entering  this  number.  Input  is  free  field : numbers 
can  be  specified  in  any  normal  way.  Should  a  mistake  be 
noticed  after  the  return  key  has  been  hit,  exit  the  program 
with  a  Control  Y  and  restart  the  program. 

Program  asks: 

INPUT  X-CENTER  OF  BUOYANCY (IN) 

Respond  by  entering  the  proper  value.  The  program  will  continue 
to  ask  for  all  the  variables  described  previously  until  after 
the  value  of  CDAT  for  the  last  component  has  been  entered.  At 
this  point  the  program  asks: 

TYPE  THE  NUMBER  3  TO  HAVE  HARD  COPY  OUTPUT 

Respond  by  typing  a  3  if  paper  tabular  output  is  desired. 

If  only  plots  are  wanted  then  type  any  other  character  and 
return.  The  program  responds  by  typing: 

INPUT  PITCH  ANGLE (DEGREES) 


After  inputting  an  appropriate  value,  the  program  will  commence 
typing  the  computed  values  of  forces  and  moments  described  above 
if  hard  copy  output  was  specified.  The  format  is: 


RUN  NO.  1 

THETA=  _ 

JDEGREES 

TERMINAL  VELOCITY 

=  _ FT/SEC 

VELOCITY 

VELOCITY 

X -FORCE 

Z -FORCE 

PITCH 

RIGHTING 

TENSION 

MOMENT 

MOMENT 

(FT/SEC) 

(M/MIN) 

(LBS) 

(LBS) 

(FT-LBS) 

(FT-LBS) 

(LBS) 

After  all  the  results  have  been  printed  out,  if  the  hard  copy  option 
was  specified,  or  directly  after  inputting  the  pitch  angle  if  no  hard  copy 
was  specified,  the  program  responds  with: 

TYPE  THE  NUMBER  1  TO  MAKE  ANOTHER  RUN 

? 

If  another  run  is  wanted  the  program  will  ask  for  a  new  pitch  angle 
and  the  above  sequence  will  be  repeated.  If  no  more  runs  are  desired 
the  program  will  respond: 

TYPE  THE  NUMBER  2  TO  HAVE  PLOTS  MADE 

? 


Respond  by  typing  a  2  if  plots  are  needed.  Any  other  character  will 
end  the  program  If  a  2  is  typed,  the  program  asks 


TYPE:  1 

2 

3 

4 


TO  PLOT  TANGENTIAL  FORCE 
TO  PLOT  NORMAL  FORCE 
TO  PLOT  PITCHING  MOMENT 
TO  PLOT  RIGHTING  MOMENT 
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Respond  by  typing  one  of  these  four  numbers  corresponding  to  the 
desired  plot.  The  program  will  then  ask: 

TYPE  PLOT  TITLE  (40  CHAR. MAX) 

Input  an  appropriate  plot  title.  At  this  point  the  program  returns  with 

TYPE  5  TO  HAVE  ANOTHER  PLOT  MADE 

Respond  by  typing  in  a  5  if  more  plots  are  wanted  or  any  other 
character  if  no  plots  are  desired.  Program  ends  when  no  more  plots 
are  desired. 

At  this  stage,  no  plot  has  been  drawn  if  plotting  was  desired.  What 
has  been  created  is  a  deferred  plotting  file  with  the  name  given  by 
the  user  in  the  SET  command  prior  to  running  the  CTDFLI  program.  To 
display  the  deferred  plotting  file  requires  the  use  of  another  program 
dependent  on  the  plotter  used.  At  WH01  there  are  three  different 
plotters,  Tektronix,  Versatec  and  Calcomp,  available  to  use.  The 
CTDFLI  program  has  been  written  primarily  for  use  with  the  Versatec 
and  this  will  be  the  only  plotter  discussed  in  the  remainder  of  this 
report.  The  Sigma  7  computer  group  can  answer  any  questions  about 
the  other  two  systems. 

The  steps  required  to  get  a  Versatec  plot  of  the  run  involve  submitting 
a  batch  job.  This  can  be  done  while  on-line  by  writing  a  small  disk 
file  with  the  appropriate  commands.  A  listing  of  the  required  commands 
is  as  follows: 

.'JOB  out,  user 
ILIMIT_(TIME,2) (CORE, 20) 

!MESSAGE  USES  VERSATEC 
:SET_F: 95 /Filename  ;  SAVE 
IPLOTV 

The  filename  in  the  SET  F:95  command  must  be  the  same  as  specified 
prior  to  running  CTDFLI.  The  SAVE  option  in  the  SET  command  prevents 
the  file  from  being  erased  after  the  plot  has  been  made.  The  above 
commands  should  be  placed  in  a  disk  file  under  a  filename  such  as 
PLOTVER.  To  have  the  Versatec  plots  made  while  on-line  requires  the 
following  command: 

BATC  H_PL0TVE  R 

This  command  instructs  the  computer  to  have  a  Versatec  Plot  drawn. 

The  computer  responds  by  typing: 

ID= _ SUBMITTED  10j.25  JUL  30,  '80 

WAITING:  TO  RUN. 


The  Versatec  plot  will  be  placed  in  the  users  bin  upon  completion. 
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RESTRICTIONS: 


STORAGE  REQUIREMENTS: 

SUBROUTINES  REQUIRED:  INPUT,  OUTPUT,  CALCUL,  IFPLOT,  PLOTFOURS  (User 

subroutines  with  main  program) 

AMITERM,  PLOTDFER  (from  acct  3  library) 


Device 


Card  Reader 
or  Terminal 

Line  Printer 
or  Terminal 


TIMING:  Fast 

PROGRAMMER:  Michael  F.  Cook 

Mary  M.  Moffett  (PLOTFOURS  SUBROUTINE) 

ORIGINATORS :  Michael  F„  Cook 

Michael  S.  Triantafyllou 

DATE:  July  1980 


Function  Special  Requirement 

Input  F: 105  DCB 

Output  F : 108  DCB 


A 


APPENDIX  "B" 

"Dynamic  Stability  Analysis 
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"Dynamic  Stability  Analysis" 

Definition  of  body-fixed  coordinate  system 
Principal  axes  (all  products  of  inertia  are  zero) 


Components 

Vector 

unit  vectors 

•*  -  fT* 

displacements 

>y 

roll  angle 

linear  velocities 

U.  (F.  u^- 

u 

pitch  angle 

CT' 

/  / 

— ^ 

angular  velocities 

h  l>  * 

n 

yaw  angle 

-JL 

hydrodynamic  forces 

x,y.  2 

F 

vector  moments 

K.H  N 

'  t 

'*1 

moments  of  inertia 


Angular  momentum 


subscript  o  -  initial  condition 

A  -  rtinnoo  In 
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Derivation  of  Equations  of  Motion 
(from  Abkowitz,  Appendix  I) 

Starting  with  Newton’s  laws  of  motion, 

r  =  JL  (m  u)  (0 

w 

and  expanding  in  terms  of  the  variables  on  the  previous  page,  six  general 
equations  of  motions  about  the  principal  axes  can  be  derived  and  the 
resulting  equations  are  shown  below. 


X  ^ -  or,  i(pf ft) 

y « ofy  "j&C1** p*)  +  ^Gr  C y ^ 

Z>  U<j  V‘)  +  x&(K-7)  (SJ 

/C  =  ly  f>  *  t  ) 

M  -  Ijj  +{4,-It)ry b  +  *»jz6(ui-uyf-irr-)-  (J) 

A/  .  Igr+C-Ty-tyff  +*H^(&^unu>f)-ys.(d  +  u>f-i>-rj(*) 


To  determine  the  forces  and  moments  above,  which  are  functions  of  geometry, 
velocity,  etc.,  the  standard  practice  is  to  expand  these  equations  in  a  Taylor 
series  expansion  about  a  suitable  initial  condition.  For  our  purposes,  we^ 
will  expand  about  steady  vertical  motion  (dc  £/0-  0;  *  ur0  *-  o  •  v0  r  u£  e  O; 


i 
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Also,  the  six  equations  of  motion  are  usually  divided  into  two  groups  for 
stability  purposes: 

1.  longitudinal  (symnetric)  stability  X f  A/ 

2.  lateral  (asymmetric)  stability  f  Kj 

In  our  case,  the  CTD  is  axisymmetric  (with  no  pinger)  so  that  only  4  equations 
need  be  solved.  C^sy j  The  four  chosen  are  x,z,/n,  *r. 

Prior  to  linearizing  equations  3,5,6  and  8,  the  following  assumptions  were  made: 

1.  CG  is  on  axis  of  symmetry  along  x  -  axis,  i.e.  a  •*  O 

2.  CTD  has  inertial  symmetry,  i.e.  T'y  -  ^  '  1*  ■  7, 

3.  Longitudinal  equations  C^t  ^  ")  assume<*  to  functions  of  only 

&  f  t C,  (jlFf  OL  U>-  )  <j  aH  other  terms  being  zero. 

4.  K  equation  assumed  to  be  function  of  ^  A;  i )  S  ^  ^ 

Using  the  above  assumptions  equations  3,5,6  and  8  become: 


X  * 

(u  +  Uj,  X^J1) 

(*> 

z  a 

'M 

(tb-  -  Uf  .  x^j ) 

0o) 

M  * 

z, 

j  ~  ^Xc.  ( •*-  *f) 

f*) 

/ 

K  = 

CnJ 

At  this  point,  these  equations  must  be  linearized  about  steady  vertical 
motion.  In  this  process  only  the  linear  terms  in  the  change  of  a  variable 
from  the  steady  position  are  kept.  As  an  example  of  how  this  is  done, 
consider  the  equation  for  X.  As  stated  in  assumption  3.,  X  is  assumed  to  be 
a  function  of  the  following: 

Since  all  the  variables  have  equilibrium  values  of  0,  except  «  in  steady 
vertical  motion  the  change  in  a  variable  due  to  a  perturbation  from  the  steady 
motion  can  be  written  as 

A  variable  =  variable 


for  all  variables  except 
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Hence  the  linear  terms  in  a  Taylor  expansion  of  X  take  the  form 


+ . +(| 

=  Stability  derivative  of 
with  respect  to 


Noting  that  /\  « 
equation  14  can  be 


is  zero  and  using  the  stability  derivative  notation  above, 
written  as 


X  =  x.  *  X- 


xe®  •* 


Similarly  for  the  left  hand  sides  of  equations  10-12. 

Now  we  note  that  the  stability  derivatives  indicate  the  change  in  force  or 
moment  due  to  a  differential  change  in  displacement,  angle,  velocity  or 
acceleration  with  all  other  variables  in  their  equilibrium  positions.  For 
simple  shapes,  values  of  these  stability  derivatives  can  be  derived.  For 
complicated  structures,  such  as  a  CTD,  values  of  these  derivatives  can  only 
be  obtained  by  model  testing. 

At  this  point  in  the  derivation,  the  right  hand  sides  of  equations  9-12 
must  be  linearized  about  steady  vertical  motion.  To  do  this  we  expand  each 
variable  to  be  the  sum  of  the  initial  condition  plus  a  perturbation  (for 
example  U  s  Substituting  these  values  into  equations  9-12, 

neglecting  higher  order  terms,  and  recognizing  that  in  steady  vertical  motion 
only  t/0  is  non-zero,  equations  9-12  become: 


X  s  shz  U 

('<■) 

Z*  /ni  (uy  >  ,  v&jJ 

op 

(/s) 

ii 

C'V 

•  .  ,  # 

U  r  Ua  /  A  CL  Ct,  =  O  • 

(j  -  Aci 

s id+fCe 


r 
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In  the  expansions  of  X,  Z,  M  and  K,  many  of  the  stability  derivatives  are 

zero.  Specifically,  X-  ,  Af  x»  .  A*?-  are  all  zero 

since  a  translation  in  x  or  z^will  not  change  the  equilibrium  condition 

of  steady  vertical  motion.  Similarly,  for  symmetric  bodies,  Abkowitz  has 

shown  thatXt©-;  Xlj  ^  Xf  j  and  X  <£>  are  z.ero  when 

evaluated  at  condition  of  vertical  motion  equilibrium  where  tr ^  ® 

Expanding  equations  16-19  with  the  above  simplifications  leaves  the  following 
4  equations  with  10  unknowns. 

Xu  U  4.  XH  AU  —  /Hf  (A. 


+  Au  *  z^uy  +  Zj  j  +  Z?  y  +  ?<S) {&-  Hf-kj)  (*) 


k\ii <^4?+  !tc*+  £  az*  ^  ^  *J, f(n) 

One  more  simplification  is  useful  at  this  point  to  get  the  equations  into 
their  simplest  forms.  First  we  choose  which  four  of  the  variables  to  be 
independent  with  the  remaining  6  dependent.  For  the  CTD  in  steady  vertical 
motion,  0 and  were  chosen  as  independent. 

Also,  it  should  be  noted,  that  all  the  above  equations  have  considered  only 
hydrodynamic  forces.  Forces  due  to  weight  and  buoyancy  must  also  be  included. 

These  equations  can  be  written 


X  m  (U/-  Bj  Ltr>  0 

(*A) 

Z  * 

(**) 

In  steady  vertical  motion  Q j  and  X  -  \A/  —  5 


as  it  should  be. 
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Substituting  (24)  and  (25)  into  equations  20  and  21  and  expanding  in  terms 
of  the  independent  variables,  the  linearized  dynamic  stability  equations 
of  motion  for  an  axissymetric  CTD  perturbed  from  steady  vertical  motion  are: 


With  proper  values  of  stability  derivatives  these  equations  could  be 
solved  and  the  dynamic  stability  investigated. 
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Calculation  of  Pendulum  Natural  Frequency 


Example:  Configuration  C.O 

with  horizontally  mounted 
CTD  instrument  (25"x7"{J) 


For  this  simple  estimate  of  the  pendulum 
frequency  the  following  assumptions  were 

1.  CTD  package  is  in  air  (for  package 
will  be  even  lower) 


All  mass  lumped  at  CG 

Angular  Equation  of  Motion  about  ^ -  axis  (for  small  angles) 

•*  Z 

JT  Q  *  where  ’T:  SHf  K  j  X  =  radius  of  gyration 

natural  frequency  =  m  %l~  ar  , /■ 

Vo*  V  X  V  A*- 


Using  assumption  2,  K=L*33  inches=2.75  ft 

and  tOM  =:  =  3. 4Z  «<*//!**. .  =  .  5A 

Calculation  of  Moment  due  to  vortex-shedding,  M. 

Assume  lowering  speed  of  3.83  ft/sec  (70  m/min) 

Flow  around  CTD  instrument  is  laminar 

with  a  shedding  frequency  of  1.3  Hz  «  8.2  rad/sec 

Using  a  lift  coefficient  of  0.5  the  maximum  lift  (side  force)  per 
unit  length  is  found  using 

^  m  *  K.  S  ^ 

for  the  rTD  instrument. 

The  overall  magnitude  of  this  force  is 

r-  4*3 (**/'*■ )  =  Ms 

The  moment  about  the  top  connection  is  therefore 

a  3  ( UJj  £ 

or 

/U  ,  3^.3  S.Z± 


J 
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The  equation  of  motion  becomes 

I  Q  =  36.3  t 

Prior  to  solving  this  equation  we  note  that  the  driving  frequency  of 
8.2  rad/sec  is  more  than  twice  as  large  as  the  natural  frequency  of 
3.4  rad/sec.  For  a  SDOF  system  with  ^  2^,  we  are  in  the 

mass-controlled  region  and  to  a  first  approximation  the  stiffness 
term  can  be  neglected.  Therefore,  we  are  left  with 

t£  *  zc.z  a-***  <s.z  d 


with 


we  have 


I  s  /*4  X 

6  =  .  34  cJ*,  e.z  t 


solving  this  for  Q  yields 

0  -  '34  (js.?)d  +  c,i: 

o^y1- 

Q*  »00>>r*4.e/  «s  .  23 


The  possibility  of  significant  vortex-shedding  thus  seems  minimal. 
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